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Abstract 

The radiative neutron capture on lithium-7 is calculated model independently using a low energy halo 
effective field theory. The cross section is expressed in terms of scattering parameters directly related to the 
^-matrix element. The cross section depends on the poorly known /7-wave effective range parameter r\ . This 
constitutes the leading order uncertainty in traditional model calculations. It is explicitly demonstrated by 
comparing with potential model calculations. A single parameter fit describes the low energy data extremely 
well and yields r\ f« — 1.47 fm^^ 
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I. INTRODUCTION 



Low energy nuclear reactions play a crucial role in Big Bang Nucleosynthesis (BBN), stel- 
lar burning and element synthesis at supernova sites [fTl-l3l. Besides placing constraints on our 
understanding of element formation, these low energy reactions play an important role in test- 
ing astrophysical models and physics beyond the Standard Model of particle physics. Often the 
key nuclear reactions occur at energies that are not directly accessible in terrestrial laboratories. 
Radiative proton capture on beryllium ^Be(p,y)^B is one of them — it is important for boron-8 
production in the sun, whose weak decay results in the high energy neutrinos that are detected at 
terrestrial laboratories looking for physics beyond the Standard Model. The relevant solar energy, 
the Gamow peak, for this reaction is around 20 keV [4]. This necessitates extrapolation to solar 
energies of known experimental capture cross sections from above around 100 keV. Theoretical 
input becomes necessary for this extrapolation. Effective field theory (EFT) is an ideal formalism 
for this as it provides a model-independent calculation with reliable error estimates. 

In an EFT, one identifies the relevant low energy degrees of freedom and constructs the most 
general interactions allowed by symmetry without modeling the short distance physics. The inter- 
actions are organized in a low momentum expansion. At a given order in the expansion, a finite 
number of interactions has to be considered and an a priori estimate of the theoretical error can 
be made. Establishing theoretical errors is crucial due to astrophysical demands ([US HI. A sys- 
tematic expansion of interactions is important because many processes involve external currents, 
and any prescription used in phenomenological models involve some uncertainty. As an exam- 
ple, the cross section for n{p,y)d at BBN energies was calculated within EFT to an accuracy of 
about 1% 0. Systematic treatment of two-body currents was necessary to achieve this level of 
precision, and it addressed a critical need ^ for nuclear theory input in astrophysics. 

While applications of EFT to systems with A < 4 nucleons is well developed, for A > 5 it is 
still in its infancy. However, some loosely bound systems, like halo nuclei open new possibilities. 
The small separation energy of the valence nucleons in halo nuclei provides a small expansion 
parameter for constructing a halo EFT [6]. The ^B nucleus with a proton weakly bound to the ^Be 
core by 0.1375 MeV is a halo system. Current extrapolation of the ^Be(/),Y)^B cross section to 
solar energies introduce errors in the 5 — 20% range |4l|7J[8]|. A model-independent EFT calcula- 
tion would be very useful to estimate the errors in the extrapolation. In addition, this would be an 
important step in developing EFT techniques for weakly-bound nuclei as has been accomplished 
in the few nucleon systems. Experiments such as those planned at the future FRIB dUl would ex- 
plore exotic nuclei near the drip lines where halo systems abound. Structures and reactions with 
halo EFT can serve as benchmark for phenomenological models of nuclei near the drip lines. 

In this paper we consider the low energy reaction ^Li(n,Y)^Li, which is a isospin mirror to 
^Be(p,Y)^B. The n-^Li system allows formulating the EFT for the nuclear interactions without 
the added complication of the Coulomb force. Traditionally ^Li(n,Y)^Li has been calculated in a 
single -particle approximation as a ^Li core plus a valence neutron interacting via a Woods-Saxon 
potential [10|. This approximation breaks down at higher energies when the internal structure of 
the ^Li core is probed, for example, near the threshold for ^Li(Y,^He)a which is about 0.5 MeV 
above the binding energy B ^ 2.03 MeV of the ^Li core. We treat the ^Li nucleus as point-like 
since we work at very low energies. Once the nuclear piece is calculated in EFT for the n-^hi 
system the Coulomb interaction in p-^Be can be incorporated systematically as have been done 
for proton fusion in EFT [fTTTl . The reaction ^Li(n,Y)^Li, besides being a check on the mirror 
^Be(p,Y)^B reaction, is important in inhomogeneous BBN. It impacts the production of carbon- 
oxygen-nitrogen in the early universe, and constrains astrophysical models [fT2l . We calculate the 
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Li(n,Y) Li reaction analytically and express the result in terms of parameters directly related to 
observables, thus quantifying the dominant theoretical uncertainty in the single particle approxi- 
mation. 



n. INTERACTION 

The relevant low energy nuclear degrees of freedom, here, are the point-like neutron, ^Li and 
^Li with spin-parity j^, | and 2+ respectively. At low energies the relevant partial waves in 
the incoming neutron-lithium state are 5-waves: ^Si, ^82 in the spectroscopic notation ^^^^Lj. The 
ground state is a 2+ state that is primarily the symmetric combination of the possible p-wave states 
^P2 and [IT3l . Conservation of parity implies that the reaction ^Li(n,Y)^Li proceeds through 
the electric dipole transition El at lowest order. 

It is known that the non-relativistic amplitude in the Z-th partial wave has the general form 

/J /j^'+^cotO/ — zp^'+^ 

where /j is the reduced mass of the n-^Li system with masses Mf^ and Mc, respectively, and 5/ is 

the partial wave phase shift. The term cot 5/ = — 1/a/ + r/p^/2-| has an analytic effective 

range expansion (ERE) for short range interactions. Scattering and bound state information in 
EFT are incorporated by matching the EFT couplings to the ERE parameters in the low energy 
expansion. 

For the initial 5-wave states, at sufficiently low momentum zJ^ ^ —i^ao, and one keeps only 
the first term in the ERE, corresponding to a single perturbative insertion of the leading EFT 
interaction. However, to describe shallow bound or virtual states that correspond to large scattering 
length ao ^ ro one has to expand around the I/ciq pole term and write 

.271 1 foP^ ^ 

This requires a non-perturbative resummation of a single interaction in EFT at leading order. Such 
a resummation extends the validity of the EFT to include the shallow state at momenta p ~ i/ao. 

In the n-^Li system, the scattering length is a^^^ = —3.63 ±0.05 fm (a^^^ = 0.87 ±0.07 fm) in 
the (^Si) spin-channel [14J. This corresponds to neutron momentum around 54 (227) MeV, or 
center of mass (CM) energy 2 (31) MeV. We are interested at the extremely low solar energies with 
momenta ;? ^ 54 MeV. Thus in EFT a single perturbative interaction in each of the ^52 and ^Si 
channels is required. The interaction can be resummed in the ^S2 channel if one wants to compare 
with data at CM energies ~ 1 MeV. 

The leading order interactions for 5-wave contain no derivatives. The spin-^ neutron and spin- 
I ^Li nucleus can be combined into the and ^52 states using the Clebsch-Gordan coefficient 
matrices Fi, Qij as N^FtC and N^QijC respectively. The vector index in Fi relates to the three 
magnetic quantum numbers in the spin 5=1 channel. The symmetric, traceless matrices Qij 
relate to the five magnetic quantum numbers in the spin 5 = 2 channel. We write the 5-wave 
interaction Lagrangian as 

=gW{NF,C)\NF,C)+g^^\NQ,jC)\NQ,jC) + --- , (3) 
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where a single momentum-independent interaction in each of the ^Si and channels was kept. 
The "■ ■ ■ " represents higher derivative terms that are suppressed at low energy. The 2x4 Clebsch- 
Gordan matrices are given as 

Fi = - ^(525,, Qij = - ^(52 [o,S, + OjS,] , (4) 



^ ^J_f-V3 1 \ =_J_f^ 01 \ .^fO I 




FIG. 1. is the ^Si, ^82 scattering amplitude. is the ^P2, ^Pj scattering amplitude. Double line is 
the ^Li propagator, single line the neutron propagator, dashed line the bare dimer propagator. 

The interaction in Eq. (|3]) produces a i'-wave amplitude shown in Fig.[T} It becomes a geometric 
series that is summed to give 

iJ^i^) (r,) - '^il. 

Lb) = -.2^^-J / (2;r)^-i^2_p2_,o+ 

where g^"^ corresponds to ^(^^ in the respective spin channels and \ is the renormalization 
scale. The loop integral L[p) is evaluated in the power divergence subtraction scheme [15] where 
divergences in both D = 4 and lower space-time dimensions are subtracted. Matching Eqs. ([2]) and 

fixes the EFT couplings as g^'^i^C) = (27r)/[//(X— \la>^)\. Introduction of the renormalization 
scale X allows for a systematic expansion of the different terms even though the final amplitude is 
independent of ?i [16J. In Ref. [1171 . initial state interactions using ERE was also considered. 

The ^Li nucleus in the final state of the reaction ^Li(n,Y)^Li is in p-wave. We will treat it as a 
shallow bound state similar to its isospin mirror nucleus. The EFT for a shallow />-wave bound 
state was formulated in Ref. [0 where it was shown that, unlike 5- wave, it requires not one but two 
non-perturbative EFT interactions. The consistent renormalization of loops is easily accomplished 
in the dimer formalism where four-fermion interactions are rewritten in terms of a spin-2 dimer 
and neutron-core interactions. The interactions in the and ^^2 state can be constructed by 
combining the matrices i^, Qij and the Galilean invariant velocity difference vector {vq — vj^\ 



(5) 
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into a />-wave state with total 7 = 2. We write the ;?-wave interaction Lagrangian as 



+ 7i:'f'- 



(6) 



where {]),■ , (tt, ,) is the dimer in the i^Pi) channel, and 



The interactions in l'^p^ are equivalent to the ones with only neutron-core short range interactions 
without a dimer field. In terms of Feynman diagrams, the four-fermion neutron-core interaction is 
replaced in the dimer formulation by a dimer exchange. Fig. [T] The non-perturbative iteration of 
the leading operators is accomplished by "dressing" the dimer propagator with nucleon-core loops. 
For a given spin-channel K = 1 CP2), 2 {^P2) the dressed dimer propagator, which is proportional 
to the elastic amplitude, reads 



iR 



I jmn 



I jmn 



(8) 



where ^ = y'—l/jpQ + /up^/M — zO+, M = Mf^ + Mc- Matching the EFT amplitudes to the p-wave 
ERE expansion determines the coupling pair (A^^^, h^^^). Again, only the first two ERE parameters 
are kept in the low energy expansion since EFT requires two operators at leading order. 



III. RADIATIVE CAPTURE 

The leading order capture cross section can be calculated via minimally coupling the photon 
by gauging the ^Li core momentum p ^ p + Zee A, where Zc = 3 is the ^Li core charge. The 
El contribution to the cross section comes from the diagrams in Fig. |2] The CM kinematics are 
defined as: p the core momentum, k the photon momentum and k p = cosG. Formally we take 
p ~ Y as the small scale where y = y/ljJE ^57.8 MeV is the ^Li binding momentum. Then at 
leading order the Mandelstam variable s ^ {Mm +Mc)^ = and |fc| = A:o ~ (jc^ + Y^) / (2^/). We 
get for the CM differential cross section 
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\k\ 




\P\ 



(i(|)(icos0 6^7(^8 \p\ SA-n^M^ 2pp 



The capture from the initial state 5*2 to the P2 final state (spin channel 2) dominates due to the 
larger initial state scattering length a'^Q^ > UqK The divergence in diagram (b) is canceled by {d). 
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(a) 




(&) 






(d) 



FIG. 2. Capture reactions ^Li(«,Y)^Li. Wavy lines represent photons. 



Summing over all polarizations and spins we get 



P2)\2 



X 



5|ZJ 



aMMn 



TlMr 



[p-i 



if 



ip 



(1+X)(1+X* 

3 



jo^sin / 2y 



p2_j_y2 \^p2_j_y2 



+ X+X* 



, (10) 



with the dimer polarization sum ££,y£*^ = Rijxy IfTSl and the wave function renormalization 

/j*^^)^|Zn| = 27r/|3y+ rj^'*!, where rj^"* is the effective range in the scattering amplitude. Z^j 
is defined as the residue at the pole in the dressed dimer propagator Dt^^pq.p) \A9l . The capture 



from the state to the state has the same exact expression as Eq. ( 10) except that and 



(2) 

r\ are replaced by the corresponding parameters in the spin channel 1. The differential cross 
section averaged over initial spin states is 



do _ 1 p^+y2l|f7lf('^2)p + |fM-('^2)|2 
JcosG ~ 327rM2 l^p 8 2 ' 



(11) 



taking the ^Li nucleus to be a symmetric combination (p/'2) + ^^2))/ "\/2 of final states. The total 
cross section o{p) is calculated with a straightforward integration over the angle 0. 

The parameters in o{p) can be determined from elastic n-^Li scattering data and ^Li binding 

(k) 

energy. However, the p-wave effective range r\ is not known accurately. This is the main the- 

(k) 

oretical uncertainty at this order. Changing the effective range r\ modifies the wave function 
renormalization factor and moves the cross section up or down by a constant factor. In tradi- 
tional potential model calculations, the parameters are determined by reproducing the ^Li binding 
energy. However, this does not constrain the effective range and other parameters of the ERE. 
For example, in a Woods-Saxon potential V{r) = — vo[l +exp(^^^^)]^^ different choices for the 

depth vo, range Rc, diffusiveness ac can be made to reproduce the known ^Li binding energy. This 
however produces different effective ranges, and constitutes an irreducible source of error in the 
theoretical calculations. 

Comparing the contributions to the capture cross section from the two spin channels analyti- 
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cally, we get 



^ ^ - 0.81, (12) 



using the same effective range ri in both spin channels. This ratio is close to the experimentally 
observed ratio [|20ll . From Eqs. ([10]), (fTTl) one can see that the total cross section at low energy is 



not independently sensitive to r\ and r^' . This is confirmed by our fit to data. 

In Fig. [3| we compare potential model calculations using Tombrello's pTO^, and Davids- 
Typel's ||71 parameters to EFT curves. At low energy the potential model results can be reproduced 
in EFT with a small variation in the effective range —0.46 fm^^ < ri < —0.3 fm^^ At higher 
energies they differ since potential models include ERE parameters beyond the scattering length 
and effective range. A fit to data from Ref. [21 J in the energy range ~ 2 — 700 eV gives an 
effective range ri = —1.83 fm^^ with only the spin channel 2 contribution and ri = —1.47 fm^^ 
with both spin channels 1 and 2. Both the ri values are compatible with the Wigner bound [22l| 
which, for a nucleon-core interaction shorter than 3 fm restricts ri to be smaller than around — 1 
fm^^. Following Ref. [|2T|. their data and the theory curves in the right panel in Fig. [s] were 
divided by the known experimental branching ratio 0.89 to the ground state and compared to a 
few other available data Il23l - l25l . The ri was fitted to the unsealed data for transition to the ground 
state as appropriate. It is clear that the theory error in the low energy extrapolation comes from 
the uncertainty in the effective range at leading order. 




FIG. 3. Potential model curves: (blue) long-dashed curve from Davids-Typel fj\, (red) dashed curve from 
Tombrello ifTOl . Left panel: (black) solid curve EFT with ri = —0.46 fm^', (black) dot-dashed curve EFT 
with ri = —0.3 fm^^ Right panel: (black) solid curve EFT with ri fitted to data. 



IV. CONCLUSIONS 

We considered radiative capture reactions for halo nuclei. The low energy ^Li(n,y)^Li cross 
section was calculated at leading order using EFT. In the single particle approximation, the cross 
section was derived in terms of scattering parameters that are directly related to 5-matrix elements. 
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Using a model-independent formalism we demonstrated and quantified the theoretical uncertainty 
associated with phenomenological potentials in the single particle approximation. The leading 
order result depends on the p-wave effective range parameter ri that is poorly known. Without 
detailed knowledge about this parameter, model calculations deviate from data at low energy. We 
extract the effective range ri by fitting our analytic form to data. 

At higher order in the EFT expansion, the cross section would get corrections from two sources: 
higher order initial and final state interactions, and two-body currents. The initial and final state 
interactions can be related to the ERE. At the very low energy, it is the final state interactions, 
which modify the wave function renormalization constants, that are important. At next-to-next- 
to-leading order the shape parameter associated with p-wave interaction contribute [l6l|26l. In 

addition, at higher order two-body currents such as Ei{NFjC)^ [NFx{V /Mc— V /Mt^)yC\Rijxy, 
where £, is the electric field, contribute. These operators are not constrained by elastic scattering. 
A higher order EFT calculation would reduce theoretical errors though at the expense of addi- 
tional parameters. This is not necessarily a drawback as what we gain is a model-independent 
understanding of the sources of higher order contributions, and a more detailed knowledge about 
the kind of experimental input that is required to better constrain the low energy theory. 

Coulomb interactions in p + ^Be scattering and ^Be(/>,Y)^B reaction is being considered where 
the current formulation plays a crucial role [26|. The power counting of electromagnetic currents 
beyond leading order is being considered as well. 
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